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ABSTRACT 

We have performed iV-body simulations of the formation of hyper- velocity stars (HVS) 
in the centre of the Milky Way due to inspiralling intermediate-mass black holes 
(IMBHs). We considered IMBHs of different masses, all starting from circular orbits 
at an initial distance of 0.1 pc. We find that the IMBHs sink to the centre of the 
Galaxy due to dynamical friction, where they deplete the central cusp of stars. Some 
of these stars become HVS and are ejected with velocities sufficiently high to escape the 
Galaxy. Since the HVS carry with them information about their origin, in particular 
in the moment of ejection, the velocity distribution and the direction in which they 
escape the Galaxy, detecting a population of HVS will provide insight in the ejection 
processes and could therefore provide indirect evidence for the existence of IMBHs. 

Our simulations show that HVS are generated in short bursts which last only 
a few Myrs until the IMBH is swallowed by the supermassive black hole (SMBH). 
HVS are ejected almost isotropically, which makes IMBH induced ejections hard to 
distinguish from ejections due to encounters of stellar binaries with a SMBH. After the 
HVS have reached the galactic halo, their escape velocities correlate with the distance 
from the Galactic centre in the sense that the fastest HVS can be found furthest away 
from the centre. The velocity distribution of HVS generated by inspiralling IMBHs is 
also nearly independent of the mass of the IMBH and can be quite distinct from one 
generated by binary encounters. Finally, our simulations show that the presence of an 
IMBH in the Galactic centre changes the stellar density distribution inside r < 0.02 pc 
into a core profile, which takes at least 100 Myrs to replenish. 

Key words: globular clusters: general - black hole physics - stellar dynamics 



1 INTRODUCTION 

Hills (1988) was the first to show that the ejection of stars 
with velocities > 1000 km s~^ is a natural consequence of 
galaxies hosting supermassive black holes. He named these 
stars 'hyper- velocity stars' (HVS). Recently, several HVS 
have been discovered in the Galactic halo (iBrown et alJ 
l2005l: iHirsch e t al."2005l: iBrown et ai]|2006allbl) . Except for 
one star which might have been ejected from the LMC 
jEdelmann et all2005ll . the travel times of all HVS are short 
enough that the stars could have been ejected from the 
Galactic centre within the lifetimes of the stars, confirming 
Hills' predictions. 

The exact formation mechanism of HVS is how- 
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ever still a matter of debate. The ejection of stars by 
super nova explosions in close bin ary systems jBlaauw 
'1961*) and dynamical encounters jPoveda. Ruiz fc Allen 
1196781 cannot produce main-sequence stars with veloc- 
ities exceeding a few hundred kilomete rs per second 
iGualandris. Portegies Zwart fc Sipiorl2o'o3) . leaving the in- 
teraction of stars in galactic nuclei around super-massive 
black holes (SMBHs) as the only possible source for HVS. 

Yu & Tremaine (2003) considered three processes which 
could eject stars from the vicinity of SMBHs: (1) close en- 
counters between two single stars, (2) encounters between 
stellar binaries and the central SMBH and (3) encounters 
between single stars and a massive black hole binary. In 
the case of the SMBH in the Galactic centre, they found 
that close encounters between single stars can eject stars 
with a rate of only 10~^^ yr~^ which would create less than 
one HVS during the lifetime of the Milky Way. The other 
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two processes were found to eject stars with rates of up to 
10"'* yr~^, sufficiently high to explain the observed number 
of HVS in the halo of the Milky Way. 

Similar results were al s o o btained by 
iGualandris. Portegies Zwart fc Sipioil (l200,'Th . who studied 
the ejection of stars by SMBHs by means of scattering 
experiments and found that HVS formation is possible by 
both processes. The tidal disruption of binaries was found 
to create HVS with higher velocities while the ejection from 
binary bla ck hole systems cr e ates H VS with a higher rate. 
Similarly, lOinsburg fc LoeTI 1)200(^1 found that the tidal 
breakup of stellar binaries can create HVS with velocities 
significantly higher than what has been found so far. 

A distinction between the two scenarios might come 
from a detailed analysis of the spatial and kinematical dis- 
tribution of HVS: the ejection of stars due to the interac- 
tion of stellar binaries with an SMBH should be nearly con- 
stant with time since the reservoir of binary stars in the 
Galactic centre is depleted only slowly. Furthermore, the 
distribution of binary orbits should be nearly isotropic in 
sufficiently relaxed nuclei, implying that the resulting HVS 
distribution will also be isotropic. In contrast, the ejection 
of stars from an SMBH-IMBH binary should show charac- 
teristic variations with time and spatial direction: HVS are 
mainly ejected when the inspiralling IMBH reaches the cen- 
tre since the density of stars is highest in the centre and the 
velocity dispersion is also highest close to the SMBH. In ad- 
dition, escaping stars acquire their extra velocities mainly in 
the direction of motion of the IMBH, introducing a spatial 
anisotropy in the HVS distribution. 

These considerations were confirmed bv iLevinI i2005l) . 
who studied analytically the distribution of escapers created 
by an SMBH-IMBH pair in a dense stellar cusp. He found 
that the ejection of stars in case of a black hole binary is hap- 
pening mainly in a burst which lasts a few dynamical fric- 
tion timescales and that if the IMBH is initially in a nearly 
circular orbit, the velocity vectors of the ejected stars also 
cluster around the orbital plane. If the IMBH moves in an 
eccentric orbit, stars are ejected in a broad jet roughly per- 
pendicular to_th«_^ainge;TeM_^^ IMBHs orbit. 
Similarly, ISesana. Haardt fc Madaul ll2006f) found through 
three-body scattering experiments that eccentric black hole 
binaries eject stars along a broad jet perpendicular to the 
semimajor axis of the binary. 

IMBHs could form in galaxies through runaway colli- 
sions of stars in star clusters, giving rise to ultra-luminous X- 
ray sources jPorteeies Zwart et al.ll2004^ . They could later 
be brought into the centres of galaxies through d^narrd - 
cal friction of the star clusters jPorteeies Zwart et al .' 2006') . 
Their merger with the central SMBHs would be an impor- 
tant source of gravitational waves, detectable with the next 
generation of gravitational wave detectors, like e.g. LISA. 
IMBHs might also be an important contributio n for the 
grow th of SMBHs in the early universe dEbisuzaki et alJ 

In the present paper we have therefore performed col- 
lisional A'^-body simulations of the dynamics of inspiralling 
IMBHs in stellar cusps around supermassive black holes. 
The aim of our simulations is to study whether the ejection 
of stars by IMBHs leads to observable consequences in the 
distribution of HVS which might help to distinguish between 
different ejection scenarios and which could give an indirect 



hint for the presence of one or more IMBHs in the centre of 
the Milky Way. 



2 DESCRIPTION OF THE A-BODY RUNS 

All runs we re performed w ith the coUisional A-body code 
NB0DY4 iAarsethI Il999l) on the GRAPE6 computers 
iMakino et alj|200ai of Bonn and Tokyo University. Our 
runs contained three different components: a central super- 
massive black hole (SMBH), an IMBH and 10^ stars. In aU 
simulations the SMBH was initially at rest at the origin and 
had a mass of Msmbh ~ 3 • 10^ M© , similar to the mass 
of the SMBH at the Galactic centre ( Schodcl ct al. 200i; 
iGhez et al.]|2005ll . The mass of the IMBHs was varied in the 
different runs. In total we performed 3 runs, using IMBH 
masses of Mimbh = 10^ Mq, 3 ■ 10^ M© and 10* M© re- 
spectively. All IMBHs moved initially in circular orbits at 
a distance of 0.1 pc from the SMBH. 

The stars had masses of m» — 30 Mq and were initially 
distributed in a cusp around the SMBHs according to the 
following density law: 

Po 



p{r) = 



^7/4 



(1) 



with Po = 3 X 10^ M© pc"^^"* and — 1 pc. For distances 
r <C ro, this distribution corresponds to a r"^/* power-law 
cusp which theoretical arguments and A-body simulations 
have shown to evolve in a stellar system around an SMBH 
(iBahcall fc Wolf .lOTg,: Baumgardt. M akino & Ebisuzak] 
l2004allbt IPreto. Merritt fc Spurzeinli2004l) . The chosen form 
for the cusp is also compatible with observations of the stel- 
lar density distribution in th e Galactic centre, w hich shows 
a power-law cusp inside 10" iGe nzel et al.ll20o3) . It has the 
additional advantage that at larger distances, where the stel- 
lar distribution has little influence on the outcome of the 
simulations, the stellar density drops off quickly. The over- 
all density was chosen in such a way to be compatible with 
current limits on th e density of stars in the Galactic centre 
jGenzel et al.ll2003ll . 

Stars were merged with the black holes (both SMBH or 
IMBH) if their separation became smaller than their tidal 
radius given by 



rt 



/Mbh\ 



1/3 



(2) 



where i?* is the radius of the stars, which we set equal to 
1 R© . The mass of disrupted stars was added to the mass of 
the black holes. IMBHs were merged with the central SMBH 
if they passed within the radius of the last stable orbit, as- 
sumed to be 3 Schwarzschild radii. We did not include the 
effects of gravitational radiation into our runs. Neither was 
any softening used in calculating the gravitational forces be- 
tween the particles in our calculation. 

Simulations were stopped when an IMBH merged with 
the SMBH or the runs reached 15 Myrs, whichever happened 
first. 

2.1 Scaling issues 

The mass of the "stars" in our simulations are 
probably too high compared to real Galactic nuclei. 
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since even if the central parts are enriched in black 
holes which s piralled into the centres through dy - 
namical f riction fBaumgardt. Makino fc Ebisxizakil (|2004h|) . 
lFreitag.~A maro-Scoanc & Kalogora (200^), the average 
mass is unlikely to be higher than a few Mq. We there- 
fore have to study to which extent the high stellar masses 
in our runs can bias our results. 

The inspiral time of the IMBHs should not change 
since dynamical friction is independent of the mass of 
the background pa,rticles as long as Mimbh ^ m, 
jBinnev fc TremaindFigSitl . which is the case both in our 
simulations and in real Galactic nuclei. Stochastic changes 
to the orbits of the IMBHs happen on a relaxation time scale 
and should scale with the average stellar mass (m) of cusp 
stars as {m)~^, i.e. in real nuclei the IMBHs change their 
orbital parameters slower than in our simulations, although 
the change should not be larger than a factor of a few. 

The number of stars removed from the cusps after the 
IMBHs have spiralled into the centre will to first order 
scale linearly with the number of stars present if the over- 
all mass density profile is fixed, i.e. would roughly be a 
factor of 10 higher in the Milky Way than in our simula- 
tions. At later stages, the inner cusps become depleted and 
stars have to be scattered into low-angular momentum orbits 
through relaxation processes in the outer cusp. The number 
of stars scattered into low-angular momentum orbits scales 
as dN ~ n{r) /tr dt, where n{r) is the number density of 
stars at radius r and tr is the relaxation time at radius r. 
Since the relaxation time scales with the mass of the stars 
as {m)~^ (Spitzer 1987), and since the number of stars n{r) 
scales as {m)~^ if the overall mass density profile is fixed, 
the number of stars scattered into low-angular momentum 
orbits is to first approximation independent of the average 
mass of stars. We will come back to this point in sec. 13.41 



3 RESULTS 

3.1 IMBH inspiral and core formation 

Initially all IMBHs sink towards the central SMBH due 
to dynamical friction. Since in all runs the masses of the 
IMBHs are much higher than the masses of the background 
stars, the frictional drag on the IMBHs is given by (see 



0.1 



iBinnev fc Tremaind lIToSTt) Eq. 7-18) 
-47r In AGV(7') A^iMBH 



dv 
'dt 



V (3) 



where p(r) is the background density of stars. In A the 
Coulomb logarithm, and X — v/{\/2a) is the ratio be- 
tween the velocity of the IMBH and the (ID) stellar ve- 
locity dispersion a. For an r~^'''^ cusp profile, a is approx- 
imately 0.60 times the circular velocity. Assuming that the 
IMBH moves in a circular orbit, and settin g InA = 6.6 
iSpinnato. Fellhauer fc Porteeies Zwartll2003ll . Eq.j^can be 
rewritten as: 



F = -47.51 



(4) 



Since the rate of angular momentum change is equal to 
dL/dt = -Fr /Mimbh and since the angular momentum itself 
is given by L = ri 
as: 



VGMsMBHf, this can be rewritten 
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Figure 1. Semi-major axis of the IMBHs as a function of time 
for the three runs performed. Solid lines show the measured semi- 
major axis, dashed lines the prediction from Eqs.|S|andQ They 
agree quite well for the initial inspiral phase. At a distance of 
about r 0.003 pc, the inspiral of the IMBHs slows down in our 
simulations because the central cusps run out of stars. 
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^SMBH 

Solving this equation with the initial condition ro = 0.1 pc 
gives for the radius reached at time t: 

4 



r{t) 



,1/4 



23.8 



Gpo Af IMBH 
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(6) 



-'SMBH 

and for the time required to reach the centre of the galaxy: 



^fric 



0.0421 



6.11 



l/4,,.3/2 
'O "-'SMBH 



^Gpo Mimbh 
MsMBH \ / Mimbh \ 



3 ■ 106 M0 y \^ 103 Mq J 



Myrs . (7) 



Fig.Qcompares the inspiral predicted by the above the- 
ory with the results of the A'^-body runs. We obtain reason- 
able agreement with the A'^-body data as long as the distance 
of the IMBHs is larger than r > 0.003 pc: the difi^erences be- 
tween the theoretical curves and the data are within the 
err ors with which the Coulomb logarithm was determined 
bv lSpinnato. Fellhauer fc Porteeies Zward i2003fl . 

Inside r = 0.003 pc the stellar cusps contain only few 
stars, (for the cusp profile chosen in our runs only M(< r) = 
2 • 10^ Mq in stars were inside this radius), so the mass in 
stars becomes comparable to the mass of the inspiralling 
IMBHs. In this case dynamical friction becomes inefficient 
and the inspiral stalls since the IMBHs cannot displace 
enough stars to lose further orbital energy. The remaining 
evolution until the simulations were stopped at T = 15 Myrs 
is mainly driven by relaxation between stars at larger radii, 
due to which some stars are scattered into the central loss 
cone around the SMBH. The resulting interactions between 



4 H. Baumgardt, A. Gualandris and S. Portegies Zwart 




5 10 
T [Myrs] 




Figure 2. Evolution of the orbital eccentricity of the inspiralling 
IMBHs. The orbits remain nearly circular during the inspiral 
phase and become highly eccentric in the stalling phase. The 
dashed lines show the eccentricity at which the inspiral time be- 
comes smaller than 10 orbital times. The IMBHs reach such high- 
eccentricity orbits within a few Myrs, implying that the lifetime 
of an IMBH in the galactic centre is very short. 



the stars and the IMBHs cause a much slower inspiral of th e 
IMBHs (Bcg elman. Blandford fc Reejll980HMakindll99'it . 

Fig-IHdepicts the evolution of orbital eccentricity of the 
IMBHs. The orbital eccentricity of an IMBH is influenced 
by dynamical friction, which arises due to many distant en- 
counters, and a few close interactions. During the inspiral 
phase, the orbits of all IMBHs stay nearly circular, despite 
a decrease in semi-major axis by almost a factor of 100. This 
is in agreement with analytic estimates which predict that 
dynamical friction in power-law halos with an isotropic ve- 
locity distribution should circular ise the orbit of an infalling 
body jTsuchiva fc Shimadallioool) . 

After the IMBHs have reached the centre and removed 
most stars from the inner cusp, dynamical friction is unim- 
portant and the orbits of the IMBHs change due to inter- 
actions with passing stars coming from larger radii. In this 
phase, the orbits acquire eccentricities as high as 1 — e = 
10~^ or larger, which, in case of the Mimbh ~ lO'* Mq 
IMBH, was high enough for the IMBH to pass within the 
radius of the last stable orbit around the SMBH, resulting 
in the merger of the two black holes. 

In the absence of perturbations, the time for two orbit- 
ing black holes to merge due to the emission of gravitational 
waves can be approximated with (Peters 1964): 



Taw = 0.0353 



G^mim2{mi + 7712) 



1 



(8) 



Here, a and e are the semi-major axis and eccentricity of 
the orbit of the two black holes, mi and are their masses 
and c is the speed of light. The dashed lines in Fig.|5|show 
the eccentricity for which the time for merging due to GW 



Figure 3. Direction of the orbital angular momentum of the 
M = 3 ■ 10^ M0 IMBH. Shown are the x (bottom panel), y 
(middle panel) and z component (top panel) of the angular mo- 
mentum vector relative to the total angular momentum. Initially, 
all IMBHs moved in the x-y plane. After the M = 3 ■ 10^ Mq 
IMBH has spiralled into the centre, the angular momentum vec- 
tor changes its direction on a timescale of a Myr or less. The 
direction of motion of ejected stars will therefore also show cor- 
relations only on small timescales while the general distribution 
should be nearly isotropic. 



emission becomes smaller than 10 orbital periods. In this 
regime the emission of gravitational waves dominates the 
evolution of the orbit, and the two black holes are likely to 
merge before dynamical interactions can reduce the eccen- 
tricity again. It can be seen that the lifetime of an IMBH in 
the galactic centre is limited to a few Myrs. The time over 
which stars are ejected by an IMBH is therefore also limited 
to a few Myrs, i.e; the majority of HVSs are generated in 
short bursts. This sets IMBH induced ejections apart from 
encounters of stellar binaries with a SMBH, which would 
create escapers nearly continuously and is therefore an im- 
portant criterion to distinguish between the two cases. 

As explained earlier, another possibility to distinguish 
the ejection of HVS due to IMBHs from binary induced 
ejections is by studying the spatial distribution of ejected 
stars. Since stars ejected by an IMBH acquire their veloc- 
ities mainly in the direction of motion of the IMBH, they 
should be ejected preferentially in the orbital plane of the 
IMBH. If the orbit is eccentric during the inspiral phase, 
stars should be ejected mainly in one direction sin ce the 
density of stars increases strongly towards the centre iLevinI 
I2OO5). In both cases it is however necessary that the orbital 
angular momentum vector and the Runge-Lenz vector of the 
IMBH stay constant for a sufficiently long time interval. 

Fig. El depicts the direction of the orbital angular mo- 
mentum vector of the Mimbh = 3 • lO'^ M© IMBH as a func- 
tion of time. Initially the IMBH orbits in the x-y plane, so L 
points towards the z-direction. For the AfiMBH ~ 3 • 10"^ M0 
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Figure 4. Density distribution of stars at three different times 
for the MiMBH = 10* Mq and Mimbh =3-10^ Mq runs. In both 
runs, the initial density distribution follows a p ~ J.—1-75 power- 
law cusp which is depleted in the central parts and turned into a 
core profile after the IMBHs have spiralled into the centre. After 
the Mimbh = 10"* Mq merged with the SMBH, the central cusp 
is replenished only very slowly (dotted lines). The density profile 
in the Galactic centre should show a similar core if it contained 
an IMBH within the last 100 Myrs. 

IMBH, the orientation of angular momentum stays approx- 
imately constant during the inspiral phase, but changes 
rapidly once the IMBH has reached the inner cusp and un- 
dergoes frequent interactions with cusp stars. After reach- 
ing the centre, the angular momentum vector changes on 
timescales of 1 Myr, and even more rapidly when the IMBH 
moves on a high-eccentricity orbit. Any correlation between 
the escape direction of different HVS is therefore destroyed 
on a timescale of a few Myrs and the general distribution 
should be nearly isotropic. We obtain this behaviour for all 
three IMBH masses and we will come back to this point 
when discussing the spatial distribution of HVS in i; l3.3l 

Fig.|l]shows the density profile of the stellar cusp for the 
Mimbh = 10* Mq and Mimbh = 3x10^ Mq runs at the start 
and by the time the runs were stopped. The initial density 
profile follows an a = 1.75 power-law cusp for radii r < 1 pc 
down to about r — 3 x 10~* pc, at which radius the cusp 
runs out of stars. Due to the merger of cusp stars with the 
central SMBH and the ejection of stars by the IMBHs, the 
initial cusps are turned into core profiles with a core radius 
of about r = 0.02 pc. The core radii stay nearly constant in 
time since the ejection rate of stars is low after the initial 
peak. The only exception is the AfiMBH ~ 10'' Mq run which 
did not last long enough to completely deplete the central 
cusp. 

Although so far no indication for a cored density pro- 
file has been found in the Galactic centre, present observa- 
tional techniques are only about now reaching the spatial 
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Figure 5. Escape rate of stars ejected from the central cusp as 
a function of time (solid and dotted lines). The stellar escape 
rate rises once the IMBHs have reached the centre and drops as 
the IMBHs deplete the central cusps. The dashed lines show the 
contribution of star-star interactions to the escape rate, which are 
responsible for only a small fraction of the escapers. Dotted lines 
show the part of the inspiral which would not have been reached 
if we had included GW emission. 

resolution to measure the density profile inside 1" . If future 
observations reveal a cored density profile this would be ad- 
ditional evidence for an IMBH in the Galactic centre. In or- 
der to test how quickly the core is replenished, we continued 
the AfiMBH = 10* Mq run after the IMBH merged with the 
SMBH. Even after running the simulation for 10 Myrs, there 
was only a slight increase in central density. Since the relax- 
ation time in our runs is about a factor of 10 smaller than 
in real Galactic nuclei, replenishing the cusp in the Galac- 
tic centre should take at least T ~ 100 Myrs and possibly 
even longer. Similar large refilling times were also obtained 
by Wang & Merritt (2005) through analytic estimates. If 
observations show that the power-law cusp in the Galactic 
centre extends down to radii much smaller than 0.02 pc, the 
presence of massive IMBHs in the Galactic centre within 
the last 100 Myrs could be excluded. In this case the ob- 
served HVS would not have been ejected by IMBHs, since 
typical travel times of HVS are of order 10* years or less 
jBrown et al.ll2006a^ . 

3.2 Ejection of HVS 

As the IMBHs sink towards the Galactic centre, interactions 
with the background stars that orbit the SMBH result in a 
large number of ejections. We consider as escaping stars all 
stars that acquire positive energies during the calculation, 
i.e. which can leave the inner cusp region modelled in our 
simulations. Fig.|S]depicts the escape rate of all ejected stars 
as a function of time. The overall evolution is very similar 
for the three IMBH masses: The escape rate rises with time 
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and reaches a maximum when the IMBHs have spiraUed into 
the centre. It drops on a similar timescale during which the 
IMBHs scatter away all stars from the inner cusp region (see 
also Fig.|lJ. At later times, escapers are created mainly by 
relaxation processes in the outer cusp, due to which stars 
with large semi-major axes can be scattered to low angular 
momentum orbits and reach the inner cusp. 

As shown in Fig.|21 the orbital evolution of the IMBH 
around the SMBH may be terminated by the emission of 
gravitational waves before the end of the calculation. We 
have therefore shown the further evolution of the escape rate 
in Fig.|^with the dotted lines. For the higher mass IMBHs 
(MiMBH ^ 3 X 10'^ M0), the majority of escapers are ejected 
well before the two black holes merge due to GW emission. 

The dashed lines in Fig.|S]show the contribution of star- 
star interactions to the escape rate. Since we did not check 
for collisions between the stars during the simulations, the 
rate of star-star escapers should be higher in our runs than 
it would be in reality. We can estimate the rate of star- 
star escapers from the Mimbh = 10'^ M0 run where in the 
beginning the IMBH moves through a low-density environ- 
ment and encounters only few stars. We find about 40 esca- 
pers/Myr in this phase, which gives an upper limit for the 
rate of star-star escapers. We find a similar escape rate for 
the Mimbh = lO" Mq run after the IMBH merged with the 
SMBH. Hence, star-star interactions contribute only a small 
fraction to the overall escape rate once the IMBHs have 
reached the centres and can therefore be neglected when 
discussing the properties of HVS. 

In order to find the HVSs in our simulations, we followed 
the orbits of all escapers in the potential of the Galaxy after 
they left the Galactic centre. We chose a fifth order Runge- 
Kutta method with adaptive step-size as the integrator. The 
model for the Galactic potential was a combination of four 
separate components: the g alactic centre, repr esented by a 
power law density profile (iGenzel et alj 1200311 . the bulge, 
represented by a Plummer model, the disc, represented by 
a Kuzmin axisymmetri c profile and the halo, represented by 
the Paczynski model (|Paczvnsdll990h . We computed the 
trajectories of all escapers for 200 Myr or until the stars 
reached a distance of 100 kpc from the centre and consid- 
ered as HVSs only those stars which acquired large enough 
velocities to be unbound and escape the Milky Way poten- 
tial. 

Fig-El shows the rate of HVS ejected from the runs. 
The escape rate of HVSs follows a trend very similar to the 
overall escape rate. It reaches a sharp maximum when the 
IMBHs reach the centre and drops as the central cusps are 
depleted in stars. 

The largest escape rates, which are reached when the 
IMBHs have spiralled down to a radius of about r — 0.002 pc 
are 150 stars/Myr for A/imbh = 10'' M©, 42 stars/Myr for 
Mimbh =3-10^ Mq and 15 stars/Myr for Mimbh = 10^ M©, 
scaling roughly with A^imbh. These numbers could be a fac- 
tor of 10 to 30 higher for real Galactic nuclei, depending 
on the central density of stars. The average rate of HVS 
after the cusps have been depleted are 15 stars/Myr for 
Mimbh = 10* Mq, 10 stars/Myr for Mimbh = 3 ■ 10^ Mq 

and 4 stars/Myr for Mimbh = 10^ M©, scaling roughly with 
1/2 

Mjj^gjj. As explained in section 2.1, we expect these num- 
bers to be independent of the average mass of stars in our 
runs. As already explained before, an IMBH creates HVS 
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Figure 6. Number of HVS created per Myr as a function of 
time for different IMBH masses. The escape rate of HVSs shows 
a similar trend to the overall escape rate and peaks when the 
IMBHs arrive in the centre. 



only in a short span of time before merging with the cen- 
tral SMBH. However, even if the IMBH can somehow avoid 
merging, HVS will show a peak in their escape times since 
depending on which IMBH mass we assume, about as many 
stars escape within 1 Myr after the IMBH has reached the 
centre as would escape within the following 40 to 100 Myrs. 

Based on the above numbers, and assuming that the 
IMBH merges with the SMBH after 5 Myrs, we can estimate 
that IMBHs with masses of 10'' Mq, 3 • 10^ Mq or 10^ Mq 
create 1700, 900 and 600 HVS respectively. If we assume a 
standard IMF, about 3% of these stars would be O or B type 
stars. Especially for the high-mass IMBHs, the numbers are 
therefore large enough to explain the 6 observed HVS. 

Fig. |7] shows the cumulative distribution of ejection ve- 
locities of HVS for the three different IMBH runs (solid line, 
dotted line, dashed line). The left panel shows the distribu- 
tion when the HVS are ejected from the cusp, while the right 
panel shows the distribution when the HVS have reached a 
distance of J? = 100 kpc from the galactic centre. The dis- 
tribution was determined from all escaping HVSs, but in- 
troducing a cut off in time due to the possible merger of 
the IMBHs has little effect. At ejection, the velocities span 
a range between ~ 800kms~^ and ~ 5000 kms^^, with the 
minimum being determined by the Galactic escape veloc- 
ity. The distribution is very similar for the three runs, im- 
plying that the velocity with which stars escape does not 
depend strongly on the IMBH mass but rather on the lo- 
cal orbital velocity at the radius where the interaction takes 
place. We compare these distributions with those generated 
by encounter s between stellar binaries and the SMBH a s 
simulated by iGualandris. Portegies Zwart fc Sipioil i2005ll . 
We consider equal mass binaries with 3 M0 stars and two 
extreme values for the initial semi-major axis; a — 0.05 AU 
(dash-dotted line) and a = 1 AU (long dash-dotted line). 
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Figure 7. Velocity distribution of HVS for three different black hole masses (solid line, dotted line, dashed line), compared with the 
distribution of escapers generated by encoun ters between stellar binaries and the SMBH (dash-dotted line and long dash-dotted line, see 
iGualandris. Portegies Zwart &: Sipiorl l2005l) V Left panel shows the velocity distribution of HVS when they are ejected from the cusp, 
right panel shows the distribution after they have reached a distance of 100 kpc from the galactic centre. The HVS distribution created 
by stellar binaries can be significantly different from the one created by IMBHs, which is nearly independent of the IMBH mass. 



The distribution for the smallest semi-major axis, which is 
the one that generates the fastest escapers, is significantly 
different from the distribution of HVS ejected by an IMBH. 
When a large sample of HVSs will be available, it will be 
possible to discriminate between the two ejection scenarios 
based on the measured velocity distribution. We notice here 
that three-dimensional velocities are needed for this kind of 
analysis. 



As expected, they all point away from the galactic centre. 
Since HVS are created in a short burst whose duration is 
much smaller than the time required to travel into the galac- 
tic halo, the escape velocities of HVS also increase in mag- 
nitude with the distance from the centre. In case HVS are 
created by an IMBH, we therefore expect a similar corre- 
lation between escape velocity and galactocentric distance 
with the fastest HVS being found at the largest distances. 



3.3 Spatial distribution of HVS 

Fig.m shows the spatial distribution of HVS for the 
MiMBH = 3000 M0 and the Mimbh = lO" M© runs 100 Myrs 
after the start of the simulation during which the orbits of 
the HVS have been followed in the potential of the Galaxy, 
assuming that the inspiral plane of the IMBH agrees with 
the plane of the galactic disc (z=0). For the Mimbh = 
3000 M0 IMBH, the distribution of HVS is nearly isotropic 
and would thus be indistinguishable from a HVS distribu- 
tion created by encounters of stellar binaries with a single 
SMBH. For the Mimbh = 10'' Mq case, there is an over- 
density of stars around the z=0 plane. This is due to the 
fact that the inspiralling IMBH had L^/y^L^ + Ll + L'i ^ 1 
during the first Myr of the inspiral, and consequently ejected 
stars preferentially in the x-y plane. For the other runs we 
could not find such overdensities and they disappear for the 
Mimbh ~ lO'' M0 IMBH if we assume that the orbital plane 
of the IMBH does not coincide with the Galactic plane. Al- 
though it is not impossible to find a non-isotropic spatial 
distribution of HVS, we conclude that this is at least un- 
likely. 

Also shown in Fig.|H|are the velocity vectors of the HVS. 



3.4 Dependence on the assumed mass of stars 

As explained in section 2.1, the average mass of stars in our 
runs is probably too high compared to real galactic nuclei. In 
order to test which influence this might have on our results, 
we repeated the Mimbh ~ SOOOMq run, using stars with 
higher masses. In the first case we used A'^ = 5 • lO** stars 
with average mass m — 60Mq and in the second case TV = 

2.5 • 10^ stars with m — I2OM0. All other parameters, like 
the density distribution of the stars and the total mass in 
stars were held constant. 

Figs. |5] and 1101 show the results for the inspiral of 
the IMBHs and the ejection of hypervelocity stars. Since 
Mimbh » m, the inspiral of the IMBHs should not de- 
pend on the mass of individual stars but only on their overall 
density and this is confirmed by Fig. |3| since all three black 
holes spiral in at more or less the same rate. The slight dif- 
ferences visible in Fig. |U] between different runs are probably 
due to statistical effects since they become important only 
in the innermost part of the cusp where only relatively few 
stars are present and do not show a trend with the stellar 
mass used. We also find that the eccentricity evolution of 
the IMBHs is independent of the stellar mass: In all 3 runs. 
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Figure 8. Spatial distribution of all HVS created in the Mimbh = 3000 Mq run (upper plots) and in the A/imbh = 10* M0 run (lower 
plots) 100 Myrs after the start of the simulations. The orbits of all escapers are followed in a realistic 3D Galactic potential to obtain 
their spatial distribution. The Galactic disc is assumed to be in the x-y plane. The distribution of HVS after 100 Myrs is isotropic in the 
x-y plane for both IMBHs, but some overdensity can be seen in the Mjmbh = 10* Mq case around z=0. Also shown are the velocity 
vectors of each star. If HVS are created by IMBHs, the fastest stars should be found at the largest distances. 



the orbits of the IMBHs stay nearly circular in the inspiral 
phase and become highly eccentric in the staUing phase. All 
IMBHs would merge with the SMBH after a few Myrs due 
to gravitational wave emission from high-eccentricity orbits. 

Fig. IIUI depicts the ejection rate of stars in the three 
runs. As discussed in sec. 2.1, we expect the maximum ejec- 
tion rate to drop linearly with the number of stars in the 
cusp and this is confirmed by our runs since we obtain max- 
imum ejection rates of N=1000, 500 and 300 stars/Myr for 
runs with iV = 10^ 5 ■ 10* and 2.5 ■ 10* stars. Fig. ITOl also 
shows that in the later phases, when the central cusp has 
been depleted in stars, the ejection rate is nearly indepen- 
dent of the number of stars which is also in very good agree- 
ment with our theoretical considerations in 2.1. We finally 
found no significant dependence of stellar escape velocity 
with the mass of stars used, and conclude therefore that the 



results of our paper should be robust against changes of the 
average stellar mass. 



4 CONCLUSIONS 

We have performed simulations of the inspiral of massive 
black holes into the centres of galaxies and of the subse- 
quent ejection of hyper-velocity stars. We found that the 
spatial distribution of HVS is nearly isotropic and would be 
difficult to distinguish from a HVS distribution created by 
interactions of stellar binaries with an SMBH if only few 
HVS were found, as is presently the case. A better indica- 
tion comes from the escape times of HVS: our simulations 
confirm that most HVS are ejected in a short burst, lasting 
only a few Myrs for typical IMBH masses, as soon as the 
IMBH reaches the galactic centre. The ejection ends when 
the IMBH merges with the central SMBH, which should 
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Figure 9. Semi-major axes as a function of time of three IMBHs 
with mass Mjm bh = 3OOOM0 moving through nuclei composed 
of stars with different stellar masses m. The resulting inspiral is 
very similar for all three IMBHs, showing that the inspiral results 
do not depend on the average mass of stars used in our runs. 



take less than 10 Myrs. Even if merging can be avoided, the 
ejection rate of HVS is a factor of 30 to 100 lower than dur- 
ing the burst majcimum. The currently observed HVS shov y 
a broad distribution of escape times jBrovyn et aljl200^ . 
which argues against ejection due to a single IMBH, but 
would still be consistent with the inspiral of several IMBHs. 
The evidence is not conclusive yet and requires more HVS 
to be found, which should be possible with future astromet- 
ric surveys like e.g. GAIA. In case HVS are ejected by an 
IMBH, we also expect a strong correlation of escape veloc- 
ity with galactocentric distance in the sense that the fastest 
HVS can be found at the largest distances. 

Another prediction from our runs is that IMBHs deplete 
the central region in stars so that an initial cusp profile is 
turned into a nearly constant density core with core radius 
r — 0.02 pc. If an IMBH was present in the Galactic centre 
within the last 100 Myrs, such a core should still be visible in 
the stellar density distribution. If on the other hand the cusp 
profile observed at larger radii continues all the way down 
to the centre, this would be evidence against the presence 
of an IMBH in the Galactic centre within the last T ~ 100 
Myrs. 
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